Efficient Classical Simulation of Optical Quantum Information Circuits^ 



Stephen D. Bartlett and Barry C. Sanders 
Department of Physics and Centre for Advanced Computing - Algorithms and Cryptography, 
Macquarie University, Sydney, NSW 2109, Australia 
(Dated: October 30, 2002) 

We identify a broad class of physical processes in an optical quantum circuit that can be efficiently 
simulated on a classical computer: this class includes unitary transformations, amplification, noise, 
and measurements. This simulatability result places powerful constraints on the capability to realize 
exponential quantum speedups as well as on inducing an optical nonlinear transformation via linear 
optics, photodetection-based measurement and classical feedforward of measurement results, optimal 
cloning, and a wide range of other processes. 
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Quantum mechanics enables information processing 
that could not be performed classically |Q; examples in- 
clude efficient factorization and secure communication. 
Optical realizations of quantum information processes 
are particularly appealing because of the robust nature of 
quantum states of light against the effects of decoherence. 
Both discrete-variable (qubit) j^, |j and continuous- 
variable (CV) H schemes offer significant potential for 
optical quantum information processing, especially if ef- 
ficient processes can be performed that are not efficient 
on any classical device. 

Advanced techniques in linear optics and squeezing, 
represented by unitary transformations of optical states, 
are known to be insufficient to implement arbitrary uni- 
tary transformations P, [§■ In order to perform uni- 
versal quantum computation ffl, optical nonlinear pro- 
cesses (such as a Kerr nonlinearity fl) have been iden- 
tified as a necessary requirement both for qubit and CV 
schemes. Nonlinear transformations are also necessary 
for other optical quantum information processes such as 
the Bell state measurements employed in quantum tele- 
portation H. Thus, the lack of a strong optical nonlinear- 
ity with low losses greatly restricts the type of quantum 
processes that can be performed in practice. Recently, 
however, nonunitary processes such as measurement have 
been identified as a means to implement nonlinear op- 
erations. Proposals for qubit-based quantum computa- 
tion [0 and CV quantum computation [Q employ photon 
counting to induce nonlinear transformations in optical 
systems. Photon counting is an important example of a 
process that can be used to achieve nonlinear transfor- 
mations via feedforward of measurement results. Such a 
nonunitary transformation appears to enable impressive 
capabilities equivalent to nonlinear transformations. 

In order to implement powerful optical quantum infor- 
mation processing that cannot be performed classically, 
it is imperative to determine what type of processes (uni- 
tary transformations, projective measurements, interac- 
tion with a reservoir, etc.) can be used to implement non- 
linear transformations and thus perform universal quan- 
tum computation. One approach is to identify classes of 



processes that can be efficiently simulated on a classical 
computer. Under the assumption that universal quan- 
tum computation is not efficiently simulatable classically, 
such processes are insufficient to implement optical non- 
linear transformations. The Gottesman-Knill (GK) the- 
orem jl], || for qubits and the CV classical simulatability 
theorem of Bartlett et al. (BSBN) || provide valuable 
tools for assessing the classical complexity of a restricted 
class of controlled unitary transformations. 

To develop a theorem that includes measurement and 
feedforward of measurement results, positive operator- 
valued measures (POVMs) are employed. In related 
work, Knill JTo[ has incorporated certain projective 
measurements into a classical simulatibility result for 
fcrmionic systems by allowing non-Hermitian Hamilto- 
nians. We place this result in a generalized and powerful 
setting: unitary transformations, POVMs and any other 
physical process can be described in the unified formalism 
of completely positive (CP) maps. We extend the defini- 
tion of the Clifford group of unitary transformations used 
in the GK and BSBN theorems to the Clifford semigroup; 
this semigroup, expressed in the language of Gaussian CP 
maps, includes a general class of unitary transformations, 
noise processes, amplifiers, and measurements with feed- 
forward. Our result is a theorem for efficient classical 
simulation of a broad class of optical quantum informa- 
tion processes (both qubit and CV), and thus a no-go 
theorem for universal optical quantum computation or 
inducing nonlinear transformations. 

Consider an optical quantum information process in- 
volving n coupled electromagnetic field modes, with a 
single mode described as a quantum harmonic oscillator. 
The two observables for the (complex) amplitudes of the 
field mode serve as canonical operators for this oscillator. 
For a system of n coupled oscillators, the 2n canonical op- 
erators {cji,pi,i = 1, . . . ,n} satisfy [qi,pj] — iMijl, with 
I the identity. We express the 2n canonical operators 
in the form of phase space coordinates, defining Zi = qi 
and z n +i — pi for i — l,...,n. These operators satisfy 
[zi, Zj] — ihY>ij, with X the skew-symmetric 2n x 2n ma- 
trix Sjj = <5i+ n j — <5i,j+n- For a state p, the means of the 
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canonical operators are defined to be & = Tr (pZi), and 
the covariance matrix is 

Itj = Tr(p(z 4 - - £,•)) - iSy . (1) 

A Gaussian state (a state whose Wigner function is 
Gaussian and thus possesses a quasiclassical descrip- 
tion) is completely characterized by its means and co- 
variance matrix Coherent states, squeezed states, 
and position- and momentum-eigenstates are all exam- 
ples of Gaussian states. 

The Clifford group C n Q is defined to be the group 
of linear transformations of the canonical operators {z{\. 
For a system of n oscillators, it is the unitary represen- 
tation of the group ISp(2n, K) consisting of phase-space 
translations plus one- and two- mode squeezing |fl^| . A 
displacement X(a) with a a real 2n- vector acts on the 
canonical operators in the Heisenberg picture as 



X(a) 



Si + aj . 



A symplectic transformation M(A) acts as 



2n 



M(A) : z j — y z\ = ^ z k a ki , 

k=l 



(2) 
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with A = (aij) a real matrix satisfying A^SA = S. A 
general element G 6 C n can be expressed as a product 
C(a,A) = X(a)M(A). 

The Clifford group consists of unitary transformations 
that map Gaussian states to Gaussian states; however, 
unitary transformations do not describe all physical pro- 
cesses. We define the Clifford semigroup, denoted /C„, 
to be the set of Gaussian CP maps jllj] on n modes: a 
Gaussian CP map takes any Gaussian state to a Gaus- 
sian state. Because Gaussian CP maps are closed under 
composition but are not necessarily invertible, they form 
a semigroup. A general element T S K, n is defined by its 
action on the canonical operators as 



T(a,A,G) 



E 



Zkaki + a.J + fit 



(4) 



where a is a real 2n- vector, G — {gij) is a In x In real 
symmetric matrix, and A is real but no longer required 
to be symplectic. Eq. (|J) includes both transformations 
(||) and (||) plus additive noise processes 1 13 described 
by quantum stochastic noise operators {fji} with expec- 
tation values equal to zero and covariance matrix 



iy^O'kiO'lj'^kl ■ (5) 
kl 



Here, pu is a Gaussian 'reservoir' state which, in order 
to define a CP map, must be chosen such that the noise 
operators satisfy the canonical quantum uncertainty re- 
lations. This condition is satisfied if the noise operators 
define a positive definite density matrix, 



G + izZ - iA^A > 0. 



(6) 



The Clifford group transformations are recovered for G — 
0. (For further details of Gaussian CP maps, see 

The action of the Clifford semigroup on the means and 
covariance matrix is given by 

T(a, A, G) : & -> £■ = £, k a ki + at 

k 

7-y -> = E a ^lkiai 3 + 9ij ■ (7) 

kl 

Because the means and covariance matrix completely de- 
fine a Gaussian state, the resulting action of the Clifford 
semigroup on Gaussian states can be easily calculated via 
this action. 

The Clifford semigroup /C n represents a broad frame- 
work to describe several important types of processes in 
a quantum optical circuit, as follows. 

Linear optics and squeezing: If G = 0, then condi- 
tion @ demands that A must be symplectic and the 
corresponding map T(a, A, 0) is C(a, A) € C n . The map 
T(a, A, 0) is thus a unitary, invertible transformation in 
the Clifford group. 

Noise, amplification, and cloning: The Clifford semi- 
group possesses an illustrative interpretation in terms of 
a linear coupling of the system's degrees of freedom to a 
reservoir. Introduction of noise into the quantum circuit 
can be incorporated using Clifford semigroup maps with 
nonzero noise term G; such noise contributions can be 
viewed as a coupling of system modes to reservoir modes 
in Gaussian states via a beamsplitter or mode coupler. 
Linear amplifiers, including phase-insensitive and phase- 
sensitive amplifiers, are also describable in this context, 
with A characterizing the amplifying term and G the 
introduction of the associated noise |Q. Linear ampli- 
fiers can be used to perform optimal cloning of Gaussian 
states 0, , describable using Gaussian CP maps . 

Measurements in the Clifford semigroup: The Clif- 
ford semigroup encompasses a broad class of measure- 
ment. Measurements as CP maps are most easily illus- 
trated using the Kraus operator sum Jl(| . A general CP 
map £ is defined by its action on an n-mode state p 
as £(p) = ~}2 k O k pO\, where {O k } are operators on the 
n-mode Hilbert space. For a measurement process, the 
operators {O k } are the elements of a POVM, and the 
index k labels a specific measurement outcome. As an 
example, let the operators O k be projections |afe)i(afe| 
onto a complete set {|«fe)i} of coherent states for the ith 
mode. The properties of coherent states ensure that the 
resulting map is Gaussian CP. In the continuous limit 
where the sum J2k ^ s replaced by the integral J d 2 a/ir, 
it describes the projective measurements corresponding 
to eight-port homodyne detection |l^] on mode i. If the 
measurement result is discarded (i.e., the reservoir, de- 
scribing the measurement apparatus, is not observed), 
this process is described by £ mapping a pure state to a 
mixture of possible measurement outcomes. 
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However, if the measurement outcome ak is recorded, 
the corresponding CP map describing the measurement is 
different, because the density matrix is updated to reflect 
our informed knowledge of the measurement outcome. In 
our example, measuring a specific result ak leads to 



P = 



\ak)i{ak\p\ctk)i{ak\ 
Tr(\a k )i(a k \p) 



(8) 



This map is also Gaussian CP. The ith mode is left in the 
known state \a k )i, and the remaining modes are collapsed 
into the corresponding multimode state, renormalized by 
the probability of the measurement outcome ak ■ 

The projections need not be onto coherent states for 
the measurement to be in the Clifford semigroup; pro- 
jections onto squeezed states are also Gaussian CP. In 
the limit of infinite squeezing, it could also be a pro- 
jection onto quadrature-phase eigenstates via homodyne 
detection 18 1. Projective measurements onto multimode 
entangled Gaussian states Jl9| (e.g., Einstein-Podolsky- 
Rosen (EPR) states) are also Gaussian CP. Furthermore, 
measurements need not be projective; Gaussian POVMs 
can be constructed by adding ancilla modes according to 
Neumark's theorem pcfl . Finally, noisy measurement can 
be described using the Clifford semigroup by composition 
with a linear noise map as described previously. 

Conditional transformations: Clifford semigroup maps 
conditioned on classical numbers or the outcome of such 
measurements can also be described using Gaussian CP 
maps. Because of the composition property of Clifford 
semigroup maps, any Gaussian CP map conditioned on 
the outcome of a Gaussian CP measurement will also be 
Gaussian CP. 

We now present the primary result of this paper, which 
is the capability to simulate processes in the Clifford 
semigroup efficiently on a classical machine. Recall that 
any Gaussian state is completely characterized by its 
means and covariance matrix. For any quantum infor- 
mation process that initiates in a Gaussian state and in- 
volves only Clifford semigroup maps, one can follow the 
evolution of the means and the covariance matrix rather 
than the quantum state itself. For a system of n coupled 
oscillators, there are 2n independent means and 2n 2 + n 
elements in the (symmetric) covariance matrix; thus, fol- 
lowing the evolution of these values requires resources 
that are polynomial in the number of coupled systems. 

Theorem: Any quantum information process that ini- 
tiates in a Gaussian state and that performs only Clifford 
semigroup maps can be efficiently simulated using a clas- 
sical computer. These maps include (i) unitary Clifford 
group transformations (displacements and squeezing), 
(ii) linear amplification (including phasc-inscnsitive and 
phase-sensitive amplification and optimal cloning), lin- 
ear loss mechanisms or additive noise, (iii) Clifford semi- 
group measurements including, but not limited to, pro- 
jective measurements in the position/momentum eigen- 
state basis or coherent/squeezed state basis, with finite 



losses, and (iv) Clifford semigroup maps conditioned on 
classical numbers or the outcomes of Clifford semigroup 
measurements (classical feedforward). 

This theorem extends beyond the simulatability results 
of GK and BSBN. Both of these previous theorems follow 
the evolution of the Pauli operators under unitary Clif- 
ford group transformations and projective measurements 
in the computational basis, allowing for the efficient sim- 
ulation of processes initiated in the computational basis. 
By instead following the means and covariance matrix of 
Gaussian states, this new theorem allows for the simu- 
lation of a much broader class of initial states, nonuni- 
tary transformations and measurements to be included 
beyond that of GK and BSBN while still incorporating 
all of their results. 

We now consider some of the key new results of this 
theorem in terms of known processes. One is that optimal 
cloning of Gaussian states cannot be used to advantage 
for any exponential quantum computational speedup. 
Also, any Clifford semigroup transformation conditioned 
on the measurement outcome of homodyne detection 
with finite losses using Gaussian states cannot be used 
to induce a nonlinear transformation, nor can a projec- 
tion onto a multimode Gaussian state be used for this 
purpose. In terms of optical implementations of quan- 
tum computing, this theorem reveals why all previous 
schemes either propose some form of optical nonlincar- 
ity §, H g § or are not efficiently scalable Ip. 

This theorem places severe constraints on the use of 
photodetection to perform nonlinear transformations in 
both qubit [[| and CV j| realizations of optical quan- 
tum computing. For a threshold photodetector |^2| with 
perfect efficiency, the POVM is given by two elements, 
corresponding to "absorption" and "non-absorption" of 
light; these elements are 



|0)(0| 



n 



\n)(n\ 



(9) 



n=l 



where \n) are Fock states of definite photon number n. 
Photon counters are effectively constructed as arrays of 
such detectors The vacuum projection describes 

the non-absorption measurement, and the corresponding 
map describing this measurement result is Gaussian CP. 
However, the absorption outcome is not. As a result, /C„ 
elements conditioned on the no-absorption outcome of a 
photodetection measurement are in /C„, whereas trans- 
formations conditioned on the absorption outcome are 
not. Note that the same result holds for finite-efficiency 
photodetectors: such detectors can be modelled as per- 
fect efficiency photodetectors with a linear loss mecha- 
nism |l^] describable (see above) using the Clifford semi- 
group. Thus, the absorption outcome of photodetection 
and the feedforward of this measurement result is a key 
resource for optical quantum information processing. 
As an example, consider the generation of single 
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photon Fock states via parametric down conversion 
(PDC) (23). The transformation corresponding to PDC is 
two- mode squeezing; thus, the production of the squeezed 
vacuum is in the Clifford semigroup. The use of a pho- 
todetector in one mode (described by Eq. (^|)) can be 
used to post-selectively create non-Gaussian states (ap- 
proximately single photon states) in the other mode con- 
ditioned on the absorption measurement outcome; such 
a process, then, is not Gaussian CP. The no-absorption 
measurement is in the Clifford semigroup, which leaves 
the other mode in a Gaussian state (the vacuum). Thus, 
the creation or use of single-photon Fock states lies out- 
side the domain of this theorem. 

Schemes for employing photon counting for linear op- 
tical quantum computation are thus constrained by two 
results of this theorem. First, linear optics gates condi- 
tioned on the non-absorption measurement of a vacuum 
cannot be used to induce a nonlinear transformation. 
The other constraint is that any nonlinear gate employing 
linear optics and photon counting must be nondetermin- 
istic; a photon counting measurement of a Gaussian state 
could possibly result in an outcome of zero photons, and 
such a result corresponds to an efficiently classically sim- 
ulatable process. (Note that nonlinear optics, in contrast, 
can be deterministic.) 

The quantum search algorithm of Grover p4| does not 
yield an exponential speedup over classical search algo- 
rithms. It is nevertheless interesting to investigate if this 
algorithm can be efficiently simulated using the methods 
presented here. We note that linear optical implemen- 
tations p5| of Grover's algorithm satisfy the conditions 
of our theorem and are thus simulatable; however, the 
resource requirements of these implementations are not 
scalable p6| ]. It is not known if a scalable optical re- 
alization would require extra resources (i.e., an optical 
nonlinearity) . 

Our theorem for efficient classical simulation provides a 
powerful tool in assessing whether a given optical process 
(such as photon counting) can enhance linear optics to 
perform nonlinear transformations or allow quantum pro- 
cesses that are exponentially faster than classical ones. 
Algorithms or circuits employing a large class of CP maps 
given by the Clifford semigroup can be efficiently simu- 
lated on a classical computer, and thus do not provide 
any sort of quantum exponential speedup p7| . Many 
quantum optics experiments can be described in terms 
of the Clifford semigroup; thus, the challenge is to de- 
velop and exploit techniques that lie outside the Clifford 
semigroup and may be used to realize powerful quantum 
information processes in an optical system. 
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